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Moments of light cone quark density, helicity, and transversity distributions are calculated in unquenched lattice 
QCD at /3 = 5.5 and /3 = 5.3 using Wilson fermions on 16® x 32 lattices. These results are combined with earlier 
calculations at /3 = 5.6 using SESAM configurations to study the continuum limit. 


1. INTRODUCTION 

Calculation of moments of light cone parton 
distributions using lattice QCD provides the only 
known way to calculate the quark and gluon 
structure of the nucleon from first principles. 
Quenched calculations of flavor non-singlet mo¬ 
ments of parton distributions yield well-known 
discrepancies with experiment [Q , so serious effort 
has been undertaken to calculate these moments 
in full QCD|,|. To augment the calculations in 
ref. ^ at /3 = 5.6, this work reports calculations 
at /3 = 5.5 and (3 = 5.3 to study the approach to 
the continuum limit using Wilson fermions. 

Moments of the quark density, helicity, and 
transversity distributions are related to the fol¬ 
lowing matrix elements of twist-2 operators: 

= 1{PS\ Q)”-' • ■ • .D^„}b|PN 
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We calculate four moments calculable using 
zero-momentum nucleon states: the spin aver¬ 
aged momentum fraction {x)u-d, the axial charge 
(l)Ati-A(i, the longitudinal spin momentum frac¬ 
tion {x)Au-Ad, and the tensor charge (l) 5 u-< 5 d- 

2. METHODOLOGY 

To facilitate direct comparison, we follow the 
same methodology as in ref. |^, which used 
SESAM configurations at /3 = 5.6 on a 16® x 32 
lattice with unimproved Wilson fermions. In par¬ 
ticular, we calculated the ratios of three-point to 
two point functions using sequential Wuppertal 
smeared sources with the nucleon projected to 
momentum zero and Dirichlet boundary condi¬ 
tions in the time direction. Renormalization con¬ 
stants to convert from lattice regularization to the 
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Table 1 

Moments of the quark distributions 
/3 K ix) (1)au—A d (^)au—A d (l)(5u—dd 

5.3 .1665 .232(6) 1.202(25) .257(9) 1.353(21) 

5.3 .1670 .253(7) 1.225(26) .276(10) 1.312(23) 

5.5 .1592 .230(3) 1.180(12) .262(4) 1.273(9) 

5.5 .1596 .238(5) 1.140(17) .265(5) 1.277(13) 

5.6 .1560 .247(9) 1.190(16) .297(6) 1.325(13) 

5.6 .1565 .242(9) 1.223(20) .298(7) 1.287(18) 

5.6 .1570 .263(13) 1.094(25) .283(9) 1.261(26) 


MS scheme at 4 GeV^ were calculated perturba- 
tively. 

At /3 = 5.6 the smearing was optimized to max¬ 
imize the overlap of the source with the nucleon 
ground state, yielding iV = 50 smearing steps, 
and the source-sink separation was chosen to be 
12 lattice steps to minimize the statistical er¬ 
rors subject to the requirement of a well-defined 
plateau in the three-point function. 

At /3 = 5.5 to keep the physical source size 
and the source-sink distance constant, we used 30 
smearing steps and a separation of 10 lattice sites. 
For K — 0.1592 we used 413 hybrid Monte Carlo 
configurations, each separated by 20 trajectories, 
and at K = 0.1596, we used 243 configurations 
separated by 10 trajectories. 

At [3 = 5.3 we used 25 smearing steps and a 
source-sink separation of 9 lattice sites. At k = 
0.1665, we used 225 configurations separated by 
15 trajectories, and at k = 0.1670 we used 240 
configurations separated by 10 trajectories. 

3. CHIRAL EXTRAPOLATION 

Since the pion cloud contributes significantly 
to many of the operators characterizing form fac¬ 
tors and moments of structure functions, it is 
reasonable to expect large, qualitative differences 
between the behavior of matrix elements in the 
heavy quark regime where current lattice calcula¬ 
tions are carried out and in the light quark regime 
pertinent to physical pions. A physical descrip¬ 
tion of the transition between these two regimes is 
given by the following interpolation formula |^[^] , 


which incorporates the leading nonanalytic be¬ 
havior in the pion mass specified by chiral per¬ 
turbation theory: 


{x'^)u-d 
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The parameter /r specifies the physical scale of 
the three quark core that serves as the source 
term for the pion cloud, and a single value of this 
parameter of the order of 500 MeV resolves the 
discrepancy in the lowest three moments {x)u-d, 
{x‘^)u-d , and {x^)u-d as well as the nucleon mag¬ 
netic moment. The best fit to the f3 = 5.6 data 
for {x)u-d is shown by the solid line in the up¬ 
per panel of Figure 1. Of particular importance 
for the present work is the fact that this inter¬ 
polation formula is nearly constant in the regime 
in which all our lattice calculations are carried 
out. Hence, even though we only have data at 
two quark masses for /3 = 5.5 and for /3 = 5.3 we 
may study the approach to the continuum limit 
by determining the best least squares fit of the 
parameter oi at each /3. Analogous interpolation 
formulae for (l)A,i-Ad, {x)Au-Ad, and {l)su-Sd 
are presented in ref. ||^ and used to fit the data 
in the lower three panels of Figure 1. Note that, 
in contrast to the spin averaged momentum frac¬ 
tion, the three spin-dependent moments are far 
less sensitive to the pion mass because of strong 
cancellation between the the nucleon and delta 
contributions in chiral perturbation theory. 


4. RESULTS 

The data are presented in Table 1. The results 
of least squares fits of the interpolation formulae 
to the data at /3 = 5.6, 5.5, and 5.3 are shown in 
Figure 1 by the solid, dashed, and dotted lines re¬ 
spectively. For each of the four operators, we see 
that the fractional shift in these lines is small as 
one changes /3. Especially since the lattice spac¬ 
ing is quite large at /3 = 5.3, these results show 
that the 0{a) corrections for unimproved Wilson 
fermions are small and benign for the operators 
considered in this work. This provides solid ev¬ 
idence that the extrapolation to the continuum 
limit is reasonably close to the results at /3 = 5.6. 
By ruling out significant a dependence, this re- 
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suit thus further substantiates the physical argu¬ 
ment that the origin of the discrepancy between 
phenomenology and naive linear extrapolation of 
lattice moments of parton distributions must lie 
in the physics of the chiral extrapolation. 
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Figure 1. Moments of quark distributions and 
interpolation formula fits at /3 = 5.6 {circles and 
solid lines), 5.5 {squares and dashed lines), and 
5.3 {triangles and dotted lines). 























